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I We present an observation of the very late thermal pulse object V4334 Sgr (Sakurai's 

. Object) with the Infrared Spectrometer (IRS) on the Spitzer Space Telescope. The 

i-Q ' emission from 5-38 /im is dominated by the still-cooling dust shell. A number of fea- 

tures are seen in absorption against the dust shell, which we attribute to HCN and 

Q ■ polyyne molecules. We use these features to determine the ^^C/^^C ratio for the ab- 

I sorbing gas to be ~ 3.2j^^^ g; this implies that, despite the H-content of the molecules, 

C/3 . the hydrocarbon-bearing gas must have originated in material produced in the very 

' late thermal pulse. We see no evidence of emission lines, despite the recently-reported 

^ I optical and radio observations that suggest the effective temperature of the stellar 

•'^ ■ remnant is rising. 

XI ■ 
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ABSTRACT 



1 INTRODUCTION 

V4334 Sgr (Sakurai's Object) is a low-mass (~ 0.5 Mq) star 
that is retracing its post-Asymptotic Giant Branch (AGB) 
evolution along the Hertzsprung-Russell diagram following 
a very late thermal pulse (VLTP) (Herwig 2001). The star 
was discovered in early 1996 (Nakano, Benetti & Duerbeck 
1996), and was found to be at the centre of a faint planetary 
nebula (PN) of diameter 40" (e.g. Pollacco 1999). It first 
appeared as an F supergiant, possibly with a hot dust shell 
(Duerbeck & Benetti 1996). 

Analysis of the optical spectrum for mid- 1996 by As- 
plund et al. (1999) revealed that Sakurai's Object had un- 
dergone a transformation to a hydrogen-deficient, carbon- 
rich star and that the elemental abundances of a range of 
species were close to those found in the R CrB stars. They 
also used the C2 (1-0) and (0-1) Swan bands at 4740A and 
5635A respectively to determine that the ^^C/^^C ratio was 
in the range 1.5-5. Pavlenko et al. (2004) used the first- 
overtone CO bands to determine that the ^^C/^^C ratio was 



4 ± 1. This is substantially less than the value in (for exam- 
ple) the Solar System C^C/^^C = 70; Mathis 2000) and 
likely arises from second stage CNO cycling following He 
burning in the VLTP (Asplund et al. 1999). 

In 1998, Sakurai's Object was completely obscured by 
an optically thick dust shell from which (as of December 
2005) it has still not emerged. Indeed analysis of the in- 
frared (IR) spectrum showed that dust emission has been 
making an increasing contribution to the IR emission since 
1997 April (Pavlenko et al. 2002) . Since the major dust event 
of 1998, the 1-5 /im spectrum is black body-like and is that 
of the ejected dust shell. Analysis of the post- 1998 IR ob- 
servations has shown that the dust is carbon, primarily in 
amorphous form (Tyne et al. 2002; Eyres et al. 2004). The 
1-5 /im spectra during the period 1999-2001 were consistent 
with massloss that increased from ~ 0.5 x 10~^ Mq yr~^ 
to ~ 1.2 X 10~^ M0 yr"^ for a gas-to-dust ratio 200, ap- 
propriate for solar abundances (note that a lower gas-to- 
dust ratio of 75 was suggested by Evans et al. (2004) for 
the hydrogen-deficient, carbon-rich wind of Sakurai's Ob- 
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Figure 1. (a): The Spitzcr IRS spectrum of Sakurai's Object in 2005. The curve is a DUSTY fit to the data; see text for details, (b): Detail 
in the spectrum in the wavelength range 12.5-16.5 /im, plotted as fx to emphasize the depression (which is indicated by horizontal 
line) around 13 //m; much of the scatter is real, arising in ro-vibrational bands of the identified molecules. 



jcct). The maximum grain radius (in a grain size distribu- 
tion n(a) da cx a^' da, with q — 3) increased by a factor ~ 3 
(Tyne et al. 2002). 

Sub-miUimetre observations (Evans et al. 2004) show 
that the flux density at 450 jim and 850 nm was still rising 
in 2003, indicating that the massloss was being maintained; 
the spectral energy distribution (SED) resembled a black- 
body at 360 K. 

Observations of [N ii]A6548, 6583 and [Oii]A7325 in the 
spectrum of Sakurai's Object in 2002 (Kerber, Pirzkal & 
Rosa 2002) hinted that the stellar temperature had increased 
from ~ 5 500 K in 1997 to > 20 000 K. The reheating seems 
to be confirmed by recent observations at 8.6 GHz with the 
VLA (Hajduk et al. 2005), which detected free- free emission 
from the newly-ionized gas. 

We describe here observations of Sakurai's Object with 
the Infrared Spectrometer (IRS) (Houck et al. 2004) on the 
Spitzer Space Telescope (Werner et al. 2003) . 

2 OBSERVATIONS 

Sakurai's Object was observed in STARE mode with the 
IRS on 2005 April 15 UT (PID 03362; AOR Key 10840320). 
Spectra were obtained with both low and high resolution 
modes, covering the spectral range of 5-38 jim, and the RED 
peak-ap array used to precisely centre the object in the IRS 
slits. We executed five 6-second ramps for a total observation 
time (including overheads) of 1042.6 s. The spectra in the 
wavelength range ~ 5.1 — 8.2 fjxn were saturated and are not 
presented here. The spectrum was extracted from the ver- 
sion 12.3 processed pipeline data product using SPICE ver- 
sion 1.1 (Spice 2005). High resolution segments were stitched 
together without any continuum normalization. 

3 DISCUSSION 
3.1 The continuum 

The dust shell, which condensed and became optically thick 
in 1998 (Duerbeck 2002) still dominates the IR SED from 5- 
38 /im (see Fig. la). Much of the scatter in the data in Fig. la 
consists of low-level features superimposed on the dust con- 
tinuum (sec Fig. lb) , and will be discussed in detail in a fu- 
ture paper. The depression in the SED around 13 yum likely 
is due to the superposition of many ro-vibrational bands of 
hydrocarbon molecules (see §3.2.1 below). 



We compare the Spitzer IRS spectrum with data 
obtained with the United Kingdom Infrared Telescope 
(UKIRT) and the James Clerk Maxwell Telescope (JCMT) 
(from Evans et al. 2004) in Fig. 2. Compared with the 2003 
SED the peak emission has clearly shifted to longer wave- 
lengths and the peak flux density has increased in 2005. The 
dust shell superflcially resembles a blackbody with temper- 
ature T ^ 200 K, signiflcantly cooler than in 2003 (Evans 
et al. 2004), when the "blackbody" temperature was 360 K. 
The flux density f\ peaks at Amax — 11.3 /xm, implying a 
temperature for the dust "photosphere" of ~ 210 K if the 
(carbon) dust is amorphous with /3 ~ 1 (where the /3-index 
is defined such that the omissivity is oc f^), or is ~ 180 K if 
the dust is graphitic (/3 ~ 2). 

The integrated flux is 1.77 x 10^^^ W m^^, yielding 
a bolomctric luminosity L, = 2 x 10^ L© for a distance 
of 1.9 kpc (Kimeswenger 2002). As the dust shell is still 
optically thick at short wavelengths, this radiation is effec- 
tively reprocessed stellar radiation and, assuming that there 
arc no 'holes' in the dust distribution (see §4.1), a measure 
of the star's bolometric luminosity. This is somewhat lower 
than the last reported in Tyne et al. (2002), who found 
L» = 4.47 X 10^ Lq . The modest decline in L« , accompanied 
by the reported rise in T, , is consistent with the evolution 
of Sakurai's Object as predicted by Hajduk et al. (2005). 

From the IRS spectrum, the "blackbody" angular di- 
ameter (Gallagher & Ney 1976) was 85 mas at the time 
of our observation, while that deduced from the 2003 
UKIRT/ JCMT data is 55 mas. In Fig. 3 we plot the black- 
body angular diameter as a function of time; earlier values 
are taken from the data of Tyne et al. (2002). The angular di- 
ameter data are consistent with linear expansion of the dust 
photosphere at the rate 0.0293 mas day~^, with the expan- 
sion starting at JD 2450763 (1997 November 10). Compar- 
ison with the visual light curve in Tyne et al. (2002) shows 
that this epoch coincides closely with the time when the 
light curve began its deep minimum (between JD 2450760 
and JD 2450850, i.e. between 1997 November 7 and 1998 
February 5). 

We fit the Spitzer dust SED using the dusty code 
(Ivezic & Elitzur 1999) with the downhill simplex routine 
wc used to fit the 1-5 /im data for Sakurai's Object (Tyne 
et al. 2002), assuming a stellar temperature of 35 000 K (Ker- 
ber, Pirzkal & Rosa 2002; Hajduk et al. 2005). We solve for 
visual optical depth (tv), graphitic carbon fraction /, maxi- 
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Table 1. Dust shell parameters for Sakurai's Object. T* is the stellar temperature assumed for the dusty fit, n is the inner radius of 
the dust shell, R* is the radius of the star, ry is the visual optical depth through the dust shell, / is the graphitic carbon fraction, amin 
and Omax are respectively the minimum and maximum grain radius, q is the exponent in the grain size distribution and Ti is the dust 
temperature at ri. 
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Figure 2. The cooling of the dust shell since 2003. Solid curve 

labelled "360 K BB" is a fit to UKIRT (squares in the range 
1-5 /im) and JCMT (squares at 450 and 850 /tm) data obtained 
in 2003 (Evans et al. 2004); data labelled "Spitzer" are reported 
here. 



mum (omax) and minimum (omin) grain size in a distribution 
n(o) da oc da, the exponent q, the temperature Ti at the 
inner radius ri of the dust shell and the ratio n /R* , where 
Jit is the stellar radius. The fit is included in Fig. 1 and the 
derived parameters arc listed in Tabic 1. For comparison, 
Table 1 also includes the corresponding parameters for the 
DUSTY fit to the 2001 September 8 UKIRT data (Tyne et al. 
2002). 

There have clearly been substantial changes in the val- 
ues of /, Ti and Omax since 2001. While the change in / im- 
plies conversion of graphitic to amorphous carbon we should 
not, as noted by Tyne et al (2002), interpret this change too 
literally as dusty computes average optical constants for a 
grain mix. The decrease in Omax (Table 1) does however sup- 
port our contention that a significant change in the charac- 
ter of both the dust mix and distribution has occurred since 
2001. Some processing of the dust has undoubtedly taken 
place, possibly a result of the grains' exposure to the hard- 
ening radiation field (Kerber, Pirzkal & Rosa 2002; Hajduk 
et al. 2005); similar changes are seen in the newly-condensed 
carbon dust around erupting novae (Gehrz et al. 1998). 

Given the observed increase in the effective tempera- 
ture of the star (Kcrbcr, Pirzkal & Rosa 2002; Hajduk ct al. 
2005), the change in Ti may imply that the inner boundary 
of the dust shell is receding from the stellar remnant, and 
that the condensation of new dust has ceased. If L* were 
constant, then n 9S T^^ implying (for no change in dust ab- 
sorptivity) a ^ 4.5-fold increase in ri since 2001. Note also 
the large increase in ri/R,, in line with the rise in n and 
the decline in Rt as the effective temperature of the star 
increases. 




JD - 2450000 

Figure 3. Variation of blackbody angular diameter with time; 
see text for details. 

Table 2. Absorption features in Sakurai's Object; identifications 
from Aoki et al. (1999) and Cernicharo et al. (2001). 
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16.085 


HCeH 
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The increase in n is strong indication that the massloss 
rate (as measured by the dust) has declined since 2001. If 
this is the case then it augurs well for the dispersal of the 
dust and the imminent reappearance of the stellar remnant 
as the dust disperses. However, in view of the way in which 
DUSTY treats optical constants in a grain mix we again cau- 
tion against over-interpretation of the results, other than to 
note the persistent decrease in T\ , and the persistent rise in 
ri, since late 1997. 

3.2 Spectral features 

3.2.1 Absorption features 

The SED is far from smooth over the entire IRS wavelength 
range. There are several obvious absorption features super- 
imposed on the dust contirmum, at 13.688 /im, 14.025 /im, 
14.135 /im, 15.079 /im, 15.918 /im, 16.085 /im (see Fig. lb). 
The only plausible identifications for these features (see Ta^- 
ble 2) are H^^CN, H^^CN, HC3N, and the carbon chain 
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molecules (polyynes) acetylene (C2H2), diacetylene (HC4H) 
and triacetylene (HCeH). 

We use the HON features to estimate the columns of 
H^^CN and H^^CN, and temperature of the absorbing ma- 
terial. For H^^CN wc take Einstein coefficients and other 
relevant data from Harris, Polyansky & Tennyson (2002), 
and partition functions from Barber, Harris & Tennyson 
(2002); the corresponding quantities for H^^CN were kindly 
provided by Dr G. Harris (private communication). The op- 
tical depth in the HCN features was dctcrniiiiod relative to 
a local continuum in the 13-17 /jm wavelength region. 

We consider an elementary model in which the absorb- 
ing gas is homogeneous and isothermal and located in front 
of the dust shell; each HCN transition is assumed to have 
intrinsic width 10 cm~^ (^ 2 x 10~^ /im) and broadened by 
a gaussian having width corresponding to the resolution of 
the IRS. The column density-temperature (iV, T) parameter 
space was searched, independently for H^^CN and H'^^CN, 
to find the best fit; we find (sec Fig. 4) 



iV = 1.18^ 



X 10^^ cm-^ T = 400 ± 100 K (H^^CN) and 



N = 5.35tlil X 10^'^ cm-^ T = 500 ± 100 K (H"CN). 

The greatest uncertainty is in the column rather than tem- 
perature; the line profile is sensitive to the latter, while the 
former is more sensitive to the placement of the contin- 
uum. We do not expect the temperatures of the H^^CN and 
H^^CN regions to differ, and the above values of T are con- 
sistent within the errors. For the purpose of discussion we 
take T = 450 K, for which = 1.42 x 10^'^ cm'^ (H^^CN) 
and N = 4.51 x 10^® cm"^ (H^^CN; see Fig. 4). This imme- 
diately gives us the ^^C/^^C ratio as ~ 3.2+j g, consistent 
with the value 4 ± 1 obtained by Pavlenko et al. (2004) from 
fitting the first overtone CO bands in the near-IR, and with 
the VLTP interpretation of Sakurai's Object. 

We note that the 1/5 band of ^^C^^CH2 is blended with 
the corresponding 1/5 band of ^^CgHg (Cernicharo et al. 
1999) and, at the resolution of the Spitzer IRS, can not be 
used to provide an alternative estimate of the ^^C/^'^C ratio. 
Similarly, we do not expect to see the effect of ^^C — > ^^C 
substitution in the polyynes in our data. However, for T = 
450 K and assuming the ^'^C/^^C ratio deduced from the 
HCN fines, the column density for C2H2 is 9.0 x 10^^ cm~^, 
where we have taken partition functions from Gamache, 
Hawkins & Rothman (1990) and spectroscopic data from 
the HITRAN database (Rothman et al. 1998). 



3.2.2 Emission features 

In view of the increased temperature of the stellar rem- 
nant (Kerber, Pirzkal & Rosa 2002; Hajduk et al. 2005) 
we expect the appearance of IR fine structure lines as 
the wind is ionized. In the case of V605 Aql, which is 
~ 100 years ahead of Sakurai's Object in its evolution, wc 
clearly see [Nell] A12.8/im (line flux 6.0 x 10"^^ W m"^) 
and [Neiii] A15.5Atm (2.6 x 10"^® W m'^) in the Spitzer 
IRS spectrum, although the dust shell is clearly much cooler 
(see Fig. 5). However there is no evidence for the presence 
of these emission lines in the Spitzer spectrum of Sakurai's 
Object, to a limit of 3 x 10"^'' W m"^ for both [Neii] and 
[Nein], that we can attribute to increased ionization of the 
ejected gas as the effective temperature of the star increases. 




Figure 4. Fit to the HCN and C2H2 features. Solid line, ob- 
served optical depth in HCN and C2H2 lines; broken line, opti- 
cal depth for column density N = 1.42 X lO" cm-^ (Ri^CN), 
4.51 X IOI6 cm-2 (H"CN) and 9.0 x lO" cm'^ (C2H2), and gas 
temperature 450 K. See text for details. 



4 DISCUSSION 

4.1 The HCN features 

The temperature of the dust "photosphere" is ~ 200 K (see 
§3.1), so we seem to have the unphysical situation that the 
temperature of the absorbing gas (~ 450 K) is greater than 
that of the dust. Cernicharo et al. (1999) note that the lAfim 
HCN feature can appear in emission if the HCN-bearing 
material is close to the star. The fact that it is in absorption 
in Sakurai's Object indicates that this material is unlikely 
to be close to the star, while the kinetic temperature of 
450 K indicates heating by ultraviolet photons in a thin 
layer of gas. For this to occur the radiation must penetrate 
the optically thick dust shell, which is possible only if the 
dust distribution has strong asymmetry. 

The HCN column is ~ 10^^ cm~^, but his should be re- 
garded as a lower limit because, as noted above, the geome- 
try of the dust and HCN is likely to be more complex than we 
have assumed here. Indeed high resolution IR spectroscopy 
(Tync et al. 2000) around the HelAl.083 line suggests that 
the distribution of matter in the outflow from Sakurai's Ob- 
ject is not spherically symmetric, as our dusty modelling 
has assumed (again underlining that the dust properties and 
distribution are not well known). Rather the dust may be 
distributed in an equatorial 'torus', to which the He and 
HCN 'wind' is orthogonal. Although the Spitzer IRS does 
not give sufficient spectral resolution to determine a doppler 
shift for HCN, in the case of the He lAl.083 line there is clear 
evidence of outflow. Tyne et al. (2000) found outflows with 
~ 600kms"\ 

However the deduced ^^C/^^C ratio should be relatively 
insensitive to the geometry. 



4.2 The H-deficient wind 

It is remarkable that the hydrogen-deficient, carbon-rich 
Sakurai's Object should display absorption lines of H- 
containing species. The presence of H-bearing molecules in 
the Spitzer IRS spectra (Fig. lb) might imply that these 
features from material ejected by Sakurai's Object prior to 
its current hydrogen-deficient state. However the "'^^C/"'^'''C 
ratio clearly demonstrates that the absorbing material must 
be a product of the VLTP. 

It is therefore of interest to consider the origin of 



36 



24 




2 



10 11 12 13 14 15 16 17 18 19 20 
X [/im] 

Figure 5. The Spitzer IRS spectrum of V605 Aql (Evans et al., 
in preparation). 

these species. In the case of IRC+10°216, the proto-typical 
carbon-rich post-AGB star, chemistry takes place in the 
photochemistry zone and polyynes are produced via the re- 
actions: 

C2H + C„H2 — > Cn+2H2 -|- H 

(Millar, Herbst & Bettens 2000). This requires the presence 
of C2H, which is formed from C2H2 by photodissociation, 
i.e. polyynes do not form before C2H2 is present. Further- 
more, Cernicharo (2004) has shown, in the context of the 

proto-planetary nebula CRL618, that carbon-rich molecules 
are formed in the presence of a strong pliotodissociating ra- 
diation field. 

It generally seems the case therefore that the chemistry 
leading to these species requires the presence of radiation. 
Such an environment might exist in the case of Sakurai's 
Object along the 'wind' axis postulated in §4.1. However 
whether the chemistry that led to the hydrocarbons we see 
in Sakurai's Object is a relatively recent phenomenon, or 
whether it is the result of earlier chemical processes whose 
products have been 'frozen in' is unclear. 

5 CONCLUSIONS 

We have reported an observation of V4334 Sgr (Sakurai's 
Object) with the Spitzer IRS. We find that the IR SED is 
dominated by the optically thick dust shell that still en- 
shrouds the stellar remnant. The dust shell "photosphere" 
has a temperature ~ 200 K, while the temperature of the 
inner dust shell is inferred to be ~ 400 K. There is evidence 
that the condensation of new dust has ceased, and that the 
character of the dust has changed since 2001. 

There are several prominent absorption features, which 
we attribute to HCN, acetylene and polyynes. The ^^C/^'^C 
ratio is ~ 3; this is substantially smaller than the Solar Sys- 
tem value (~ 70) and is consistent with (a) the CNO cycling 
following the VLTP, and (b) the value obtained from the IR 
first overtone CO bands. 

We see no evidence for the presence of emission lines 
(such as [Neil] A12.8/im, [Neiii] A15.5/im) that would be a 
signature of the anticipated increased T^g of the star. 
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